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ABSTRACT 


A  general  formulation  of  the  boundary-layer  computation  scheme  on  the  surface  of  a 
rotating  blade  is  presented.  Momentum-integral  methods,  together  with  the  three-dimensional 
entrainment  equation  for  a  rotating  disk,  arc  used  to  calculate  the  three-dimensional  turbulent 
boundary  layer  in  an  orthogonal  streamline  coordinate  system.  First -order  finite-difference 
methods  arc  used  to  solve  the  resulting  boundary-layer  equations.  The  unknown  variables  are 
the  st Teamwise  momentum  thickness,  the  shape  parameter,  and  the  streamline  slope  at  the  sur¬ 
face.  The  boundary-layer  calculation  method  is  combined  with  existing  geometrical  and  imiscid 
potential-flow  computer  codes  for  rotating  blades  to  form  an  efficient  turbulent  boundary-layer 
computer  code.  For  a  given  potential-flow  solution,  a  typical  boundary-layer  computation  re¬ 
quires  less  than  10  seconds  computer  time  on  the  Burroughs  77()()  high-speed  computer. 
Boundary-layer  predictions  are  presented  for  several  rotating  blades.  Computed  results  are 
shown  to  be  in  agreement  wit!:  experimental  data  for  a  simple  rotating  body,  l  or  the  examples 
considered,  displacement  of  the  mid-chord  point  of  a  blade  from  a  straight  radial  line  is 
predicted  to  reduce  the  computed  values  of  local  skin  friction  coefficient,  with  an  estimated  it: 
crease  in  overall  efficiency  of  about  one  percentage  point. 
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INTRODUCTION 


Turbulent  boundary  layers  occur  in  practice  on  the  rotating  blades  ol  compressors,  turbines, 
helicopters,  or  propulsors.  Centrifugal  and  Coriolis  forces  are  additional  terms  that  arise  in  the 
governing  equations  due  to  the  blade  rotation  relative  to  the  description  of  non-rotating  systems. 
These  forces,  along  with  the  pressure  and  viscous  forces,  result  in  a  complex  three-dimensional 
How  configuration. 

Several  previous  studies  have  been  made  of  turbulent  boundary  layers  on  rotating  bodies. 
One  of  the  earliest  studies  was  by  von  Kantian*11  in  046  on  axisymmetrie  bodies.  Banks  and 
Ciadd*-1,  in  062,  investigated  both  laminar  and  turbulent  boundary  layers  on  rotating  blades.  \ 
review  of  much  of  the  recent  work  is  given  by  Meyne*-*1  (072). 

In  069 .  Cham  and  Head*-11  presented  experimental  and  theoretical  turbulent  boundary -lave 
results  on  a  rotating  disk .  In  following  years,  they  extended  their  theory  to  include  the  flows 
over  a  rotating  circular  cylinder  it  an  axial  stream*'1  and  a  rotating  axisymmetric  body**’1.  Ali 
three  of  Cham  and  Head's  investigations  used  the  momentum  integral*  1  with  entrainment*1'1 
method  of  solution,  f  he  assumed  entrainment  function  and  the  exact  form  of  the  momentum  in 
tegrai  equations  varied  with  the  differing  geometries.  Cood  agreement  was  lound  between 
theoretical  piedictioits  and  experimental  results  for  Cham  and  Head's  investigations. 

I  akshminaray  ana.  el  al . • presented  results  ol  an  in'  estimation  of  the  turbulent  houndaiy 
lay ci  on  a  rotating  helical  blade  in  l'C2.  Again,  the  solution  was  based  on  the  momentum  m 
tegrai  method.  I  lie  governing  equations  were  simplified  by  the  assumptions  ol  zero  pressiue 
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gradient  and  the  use  of  the  orthogonal  streamline  eoordinate  system.  Overall,  their  theoretical 
predictions  were  in  acceptable  agreement  with  their  experimental  results:  measured  values  of  the 
streamwise  momentum  thickness  were  no  more  than  20  percent  less  than  predictions  at 
midradius  and  values  of  the  limiting  streamline  were  in  a  2  degrees  agreement  with  predictions. 

Also  in  1972,  Meynel')  presented  theoretical  and  experimental  results  for  boundary-layer 
flow  on  rotating  wedged-slutped  surfaces.  Meyne  sought  a  generalized  similarity  solution  (reduc¬ 
ing  the  equations  to  an  ordinary  differential  equation  in  the  normal  direction)  and  considered 
both  laminar  and  urbulent  flow.  Detailed  results  are  shown  for  the  predicted  boundary-layer 
parameters  for  laminar  flow.  However,  his  equations  for  the  turbulent  flow  case  were  identical 
to  those  of  Banks  and  (.iadd*-*  and  tints  no  new  information  was  obtained. 

1  he  present  work  extends  that  of  the  previous  investigators  by  providing  lot  theoretical 
computations  of  the  turbulent  boundary  layer  that  develops  about  more  complex  geometries.  In 
particular,  the  current  method  was  designed  to  compute  the  turbulent  boundary  layer  on  the  sur¬ 
faces  of  rotating  blades  with  thickness  and  loading  in  an  axial  stream  operating  at  0.7-radius 
Reynolds  numbers  greater  than  about  4  x  1()6  lot  which  fully  turbulent  flow  is  expected. 
Although  not  considered  explicitly,  the  present  calculation  method  could  be  modified  slightly  to 
perform  the  laminar  calculation  as  well. 

I  he  thtee-dimensional  turbulent  boundary  layer  in  streamline  coordinates  U.  c)  is  govern 
ed  by  three  equations,  fliese  are  the  4  and  q  momentum  integral  equations  and  the  continuity 
equation,  fliese  equations  are  the  same  as  those  given  by  Smith"1*)  mid  by  f  'umpsty  and 
Head* 1 "  with  one  additional  term  (to  account  lot  the  rotation)  in  both  of  the  momentum  in¬ 
tegral  equations.  In  the  integral  method,  the  equations  of  motion  are  integrated  across  the  boun¬ 
dary  layer,  producing  first-order  partial  differential  equations.  I  mpirieal  information  is  incot- 
pointed  in  the  equations  to  account  for  Reynolds  stresses.  In  the  numerical  method,  cl  it  ect  use  is 
made  of  the  streamwise  momentum  thickness  B||  and  the  shape  parameter  H.  the  two 
parameters  which  most  readily  characterize  the  behavior  of  turbulent  boundary  lay ct  s(  1 2 )  \  |K- 
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I 

use  of  streamline  coordinates  simplifies  both  the  boundary  conditions  and  the  governing  equa- 

< 

tions  but  requires  the  definition  of  a  new  coordinate  system  based  on  the  How  field  and  blade 

geometry.  , 

Several  assumptions,  which  are  discussed  in  a  later  section  of  this  paper,  are  made  concern¬ 
ing  the  behavior  of  the  flow.  These  assumptions  reduce  the  number  of  unknown  quantities  in 
the  three  governing  equations  to  three  variables.  These  are  the  streamwise  momentum  thickness 
©I  ],  the  shape  parameter  H,  and  the  tangent  of  ft,  the  angle  between  the  direction  of  the  wall 
friction  vector  (or  wall  streamline)  and  the  inviseid  streamline. 

When  computing  the  turbulent  boundary  layer  on  a  rotating  blade  with  thickness  and  lilt,  a 
potential-flow  computer  code*1-1)  developed  at  DTNSRDC  by  Brocket!  is  used.  This  code  com¬ 
putes  the  potential  flow  velocities,  slopes  of  the  streamlines  and  crosswise-normal  lines,  metric 
coefficients,  and  the  angular  velocity  component  normal  to  the  surface.  All  values  are  given  in 
non-dimensional  coordinates  (xc,  xK)  where  \t.  is  the  chord  fraction,  measured  from  the  leading 
edge  and  xR  is  the  fraction  of  radius,  measured  from  the  axis  of  rotation. 

Streamlines  and  crosswise-normal  lines  are  determined  from  the  given  slopes  using  an  lulcr- 
predictor/trape/.oidal-corrector  method'1^.  Potential-flow  velocities,  metric  coefficients,  etc.  arc 
interpolated  in  the  streamline  coordinate  system  using  a  spline-under-tension  method' |S).  A  first- 
order  finite-difference  approximation' is  used  to  solve  the  system  for  the  three  basic  unknown 
quantities. 

No  interaction  with  the  hub  is  considered  in  the  present  analysis  and  the  nature  of  the  flow 
near  the  tip  is  assumed  to  be  adequately  predicted  by  the  potential-flow  model. 

Results  of  the  turbulent  boundary-layer  calculation  arc  presented  for  several  rotating 
geometries.  The  first  body  studied  is  the  rotating,  weakly  helical,  segment  investigated  ex¬ 
perimentally  and  analytically  by  1  akshminarayana.  el  al.'9'.  Computed  results  are  in  agreement 
with  the  experimental  data.  The  remaining  geometries  are  described  analytically  with  no  previous  1 

experimental  or  analytical  boundary -layer  data  available.  I  he  effect  of  blade  warp  aiul  skew  on 
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the  predicted  values  of  the  local  skin-friction  coefficient  is  investigated  qualitatively  using  these 
surfaces.  Blade  warp  is  defined  as  the  angular  displacement  of  the  midchord  point  from  the 
reference  line  in  the  plane  of  rotation  and  blade  skew  is  defined  as  the  displacement  of  the  mid¬ 
chord  point  along  the  pitch  helix  at  constant  radius. 

COORDINATE  REFERENCE  FRAMES 

rite  geometry  and  boundary-layer  characteristics  of  rotating  blades  ate  described  by  two 
coordinate  systems  fixed  in  the  rotating  surface.  Brockett* gives  the  geometry  and  potential 
flow  characteristics  in  terms  of  the  non-dimensional  coordinates  \c  and  \^.  sec  Figure  1.  1  liesc 
are  noil-orthogonal,  helical  coordinates  where  \c  is  the  fraction  of  chord  as  measured  from  the 
leading  edge  and  is  the  fraction  of  radius  as  measured  from  the  axis  of  rotation. 


Figure  1  Blade  Surface  Ceometix 


The  turbulent  boundary-layer  scheme  uses  an  orthogonal,  curvilinear  coordinate  system  bas¬ 
ed  on  the  potential-flow  streamlines,  but  with  explicit  dependence  on  the  xc  and  xR  coordinates. 
The  use  of  the  orthogonal  coordinates  simplifies  the  boundary  conditions  and,  hence,  the 
boundary-layer  equations  governing  the  turbulent  How.  The  orthogonal  streamline  coordinates 
on  the  surface  are  4,  which  varies  along  streamlines,  and  q,  which  is  normal  to  the  streamline. 
The  4  coordinate  measures  distance  along  the  surface  normal. 

In  the  orthogonal,  curvilinear  coordinates,  the  distance  between  lines  of  constant  4  or  q 
varies  from  point  to  point  on  the  surface.  Metric  coefficients,  denoted  by  h  j ,  hi,  and  I13  for  the 
4.  q,  and  4  coordinates,  respectively,  are  required  to  correlate  surface  distances  to  increments  of 
the  streamline  coordinates.  Since  the  boundary-layer  thickness  is  assumed  to  be  small  and  the  4 
coordinate  is  assumed  to  measure  distance  along  the  surface  normal,  the  metric  ld  can  be  set 
equal  to  unity  with  no  loss  of  generality.  A  general  element  of  length  dx  in  streamline  coor¬ 
dinates  is  then  given  by 

dx-  -  h,  2<J42  +  h,  2dr)2  +  ^2  (I, 

where  hj  and  hs  are  functions  of  4  and  q  only. 

BASIC  EQUATIONS 
Momentum  liquations 

The  continuity  and  momentum  equations  for  laminar  How  in  a  rotating  system  are  readily 
available  in  the  literature,  e.g.,  Moore*17*.  The  turbulent  boundary-layer  equations  are  the 
laminar  equations  plus  additional  terms  due  to  the  Reynolds  stresses.  In  streamline  coordinates, 
the  equations  for  turbulent  flow  are 

3(hiu)  3(h|V)  3(hihsw) 

C  ontinuitv:  — rr—  *  — r -  4  - r- —  "  0  (2a) 

34  3 q  34 
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h,  34  h2  3q  34  1 
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^-Momentum: 
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where  u,  v,  and  w  are  the  time-mean  velocity  components  in  the  4,  >7  and  4  directions,  respectively, 
K[  is  the  geodesic  curvature  of  lines  of  constant  4, 


L  £Jn 

h|h2 


Kl  =  —  — 1 


k->  is  the  geodesic  curvature  of  lines  of  constant  rj. 


*2 


1  3h2 

h)h2  34 


U  is  the  speed  along  the  streamline  at  the  outer  edge  of  the  boundary  layer, 

H  is  the  viscosity  of  the  fluid, 
q  is  the  density  of  the  fluid, 

cut  is  the  angular  velocity  component  in  the  4 -direction  (11-direction) 

co 3  =  co  ■  n  =  co(j_  ■  11) 

u'w'  and  v'w'  are  the  time  averaged  products  of  the  turbulent  fluctuating  velocity  components 
(i.e.,  components  of  the  Reynolds  stresses). 

The  integral  boundary-layer  equations  are  obtained  by  integrating  the  momentum  equations 
with  respect  to  4  from  the  surface  of  the  blade  (4  =  0)  to  the  edge  of  the  boundary-layer  (4=3). 
The  continuity  equation  is  used  to  eliminate  all  terms  involving  the  velocity  component  .  nor¬ 
mal  to  the  blade  surface.  The  resulting  equations*  ,l»  are: 


4-momentum  integral: 


1  30ii  1  30i2 

—  -  + - —  +  K2(0|  1  -  022)  +  K  | (20s |  -  32) 

h|  34  h2  3 17 

2021  -d2  3U  20) )  +3,  3U  2w3  c 

Uh2  3 17  Uh)  34  U  2  1 


(3a) 
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^-momentum  integral:  —  ^21  +  ~  ^21  +  2871  (  — ~  —  +  K2 

h,  34  h2  dr)  21  \Uh,  34  2 


(3b) 


2cl>3 


2822  3U 


— +T^T,  -M«l.-«22^.)-7Cf2 


In  equations  3a  and  3b,  0jj,  0]2,  02j,  022,  dj,  and  d2  are  the  momentum  and  displacement 
thicknesses  defined,  respectively,  by  succeeding  equations  5a  through  51  and 

3v 


cr.  = 


3u  — ; — ; 

^  34  ~  6U  V' 

yeu2 


and  Cjc,  = 


u -  -  £>v  w 

H  34 


yeu2 


(4) 


arc  the  turbulent  wall  skin  friction  coefficient  components  in  the  4  and  rj  directions,  respectively. 
The  momentum  and  displacement  thicknesses  are  defined  by: 


»u  -  /(■  -c)0« 


(5a) 


6)2  “  {(’  -uv)hA' 


(5b) 


uv 


02i  =  ~/-yd4 

0  L 

ds  •  8 | 1 


<5c) 


6  \2 


022  =  -y  d4 

0 


(5d) 


3)  =  /  (1  -  —  )d4 


(5c) 


d2  =  -  /  d4 
0^ 


(5f) 


The  boundary  conditions  are  u  =  v  =  0  at  4  =-=  0  and  u  U,  v  0  at  4  =  d.  the  edge  ol 
the  boundary-layer. 
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The  two  momentum  integral  equations  are  insufficient  to  define  the  boundary-layer 
development.  An  additional  equation  (the  entrainment  equation)  will  be  introduced  in  the  next 
section  and  several  assumptions  (for  Cf  ,  Cf^,  and  mean  velocity  profile  shapes)  applied  in  order 
to  reduce  the  number  of  unknown  quantities  to  three.  These  three  quantities  are  the  streamwise 
momentum  thickness,  0,i;  the  shape  factor  parameter,  H;  and  the  tangent  of  the  angle  between 
surface  streamline  and  external  streamline  directions,  tan  ft. 

Entrainment  Equation 

The  three-dimensional  entrainment  equation  of  Myring<8),  which  has  also  been  used  by 
Cumptsy  and  HeadUD  and  von  Kerczek  and  Langan<l8>,  is  the  additional  equation  used  here  to 
describe  the  turbulent  boundary-layer  flow  over  the  blade  surface.  This  equation  describes  the 
growth  of  the  boundary  layer  with  time  as  the  turbulent  mixing  process  spreads  the  general  mo¬ 
tion  of  turbulent  flowU^),  The  continuity  equation  is  integrated  with  respect  to  4  through  the 
boundary  layer  to  obtain  (in  streamline  coordinates)  the  entrainment  equation: 


2  ajcMo  .  j.  (  JL  +  K. )  (d-d,) 

h,  34  h 2  3n  '  Uh,  34  “/  1 


.  /  1  3U 

2  h2U  dr\  + 


i  a  6 


w(d) 

u 


(6a) 


34  h2  on  >  unj 

K|  )  hi  34 

The  right-hand  side  of  equation  (6a)  is  a  measure  of  the  rate  at  which  fluid  enters  the  boundan 
layer.  Denoting  this  quantity  by  E,  where  E  is  a  dimensionless  entrainment  function,  the  entrain 
ment  equation  becomes: 


i  a  (d-d,)  i  ad->  ,  i  au  .,  >  .  . . 
h;  -ir-  -  hi  ^T+  ( Uhi  ar +  K-) 


,  1  au  ..  , 

'  hsU  dri  1  ' 


h2U  3  rj 

Details  of  the  entrainment  function  are  given  in  a  later  section  of  this  paper. 


(6b) 


ADDITIONAL  RELATIONSHIPS 

f  urther  assumptions  are  needed  in  order  to  solve  the  two  momentum  integral  equations  (,Lt 
and  3b)  and  the  entrainment  equation  (6a). 
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Skin  Friction 


The  turbulent  skin  friction  terms  and  Cf2  on  the  right-hand  side  of  the  |  and  17  momen¬ 
tum  integral  equations,  respectively,  are  obtained  from  empirical  data  for  the  skin  friction  coef¬ 
ficient  Cf  in  the  local  inviscid  streamline  direction!18)  by 


Cf,  =  cf 

(7a) 

and 

Cf2  =  Cf  tan  ft 

(7b) 

where 

d\ 

V  34 

tan  p  =  hm  —  -  ltm  3u 

4-o  i—o  ^ 

(7c) 

By  assuming  that  the  component  of  boundary-layer  flow  in  the  inviscid  streamline  direction 
is  approximately  two-dimensional,  the  skin  friction  coefficient  C,  can  be  specified  using  a  two¬ 


dimensional  skin  friction  formula.  The  skin  friction  relationship  given  by  Thompson*211)  anc,  ap. 
proximated  by  Head  and  Patel*211  is  used  in  the  present  calculation  method.  The  skin  friction 


coefficient  is  approximated  by  the  analytic  relation 

Cf  =  exp  (aH  +  b)  (8a) 

where  a  =  0.019521  -  0.386768c  +  0.028345c2  -  0.00070 lc3  (8b) 

b  =  0.191511  -  0.834891c  +  0.062588c2  -  0.001953c3  (8c) 

c  In  Rgj  1  (8d) 


1)0,1 

ReM  —  (St) 


U  =  speed  at  the  edge  of  boundary  layer  (8g) 

v  =  kinematic  viscosity  of  water  ( 8 h ) 

0 1 1  -  the  dimensional  streamline  component  of  the  momentum  thickness.  (8i) 

Both  coefficients  a  and  b  decrease  as  0 j  |  increases.  Thus,  for  decreases  or  “small"  increases  in 
the  shape  factor  H,  the  skin  friction  coefficient  C,  decreases  as  0 j  t  increases. 
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Velocity  Profiles 


Assumptions  involving  the  formulation  of  the  streamwise  and  crosswise  velocity  profiles  are 
made  in  order  to  relate  the  crosswise  momentum  and  displacement  thicknesses  e12>  021’  022-  and 
d2  to  the  streamwise  momentum  and  displacement  thicknesses  0u  and  dj.  Extreme  accuracy  is 
not  necessary  in  the  representation  of  these  velocity  profiles  since  they  are  used  only  to  relate  in¬ 
tegrated  values  of  the  crosswise  and  streamwise  momentum  and  displacement  thicknesses. 

The  approximation  for  the  streamwise  velocity  profile  is  given  by  Cooked)  and  by 
Patel,--,).  This  formulation,  which  approximates  the  experimental  swept-wing  boundary-layer 
data  of  Wallace*-4*,  is 


H-l 


(9a) 


The  crosswise  velocity  profiles  are  assumed  to  be  related  to  the  streamwise  velocity  profiles 
by  empirical  formulation  due  to  Mager*--1'*.  This  approximation,  also  assumed  by  C'ooke*--*, 
Cham  and  Head*4*,  and  von  kerc/ek  and  l.angan*****,  is  given  by  the  relation 


U 


—  (l-i 

V  '  d 


tan/? 


(9b) 


This  approximation  is  assumed  to  be  adequate  as  long  as  it  is  not  applied  to  crossflow  reversal. 

An  approximation  tor  the  boundary-layer  thickness  d  is  found  by  substituting  the  assumed 
form  of  the  streamwise  velocity  profile,  equation  (9a),  into  the  definition  of  the  streamwise 
momentum  thickness  0||.  equation  (5a).  Integration  and  simplification  give  the  approximation 


J  =  (fT7  +  H) 

Crosswise  1  hick  nesses 

The  expressions  lot  streamwise  velocity  profile,  equation  (9a),  and  the  crosswise  velocity 
profile,  equation  (9b)  are  substituted  into  the  definitions  of  the  streamwise  displacement 
thickness  and  the  crosswise  momentum  and  displacement  thicknesses,  equations  (5b)  through 


(50.  and  integrated  across  the  boundary  layer.  After  considerable  manipulations,  the  following 
expressions  result: 


0 12  —  0,,  M)2  (10a) 

02i  =0n  M2i  (10b) 

022  =  0,,  M22  (10c) 

6,  =0n  D,  (100) 

d2  =0n  D2  (10c) 


whereMl2  =  -(G  +  H)tan  ft  (  +  hT2  +T TTl  *  hT?  ) 
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M21  =-(G  +  H)  tan/3  ^  ‘  H~fT  +  H  +  2  ^ 

M22  =  -(G  +  H)  tan2/3  (  —  -  H  +  ,  +  H  +  2  '  h  +  3  +  H  +4  ) 
D,  =  H 

°2  =  -(G  +  H>  tanP  (  h  +  1  '  H  +  3  +  H  +  5  ) 


(11a) 

(lib) 

(11c) 

(lid) 

(lie) 


Lntrainment  Function 

The  entrainment  function  F  is  assumed  to  be  determined  uniquely  by  the  circumferential 
velocity  profile  in  the  same  way  as  proposed  by  Head*'9'  for  a  two-dimensional  boundary  layer. 
However,  C  ham  and  Head*4*  found  that  if  the  entrainment  function  is  reduced  to  approximately 
two-thirds  of  the  Hat  plate  value,  then  agreement  between  calculated  turbulent  boundary -layer 
development  on  a  rotating  disk  and  experiments  was  excellent.  Incorporating  the  results  ot 
Cham  and  Head  into  the  entrainment  function  presented  by  Green,  et  al.,2f,>  gives  the  following 
universal  expression  for  the  function  F: 

F  =  0.025  (  -  I  )  (l2a) 
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Ci 


2H 

H-l 


(12b) 


=  the  shape  factor  parameter,  H  = 


H  =the  shape  factor  parameter,  H  =  (12c) 

The  relationship  between  the  boundary-layer  thickness  d  and  the  streamwise  momentum 
thickness  0||  given  in  equation  (9c)  can  be  rewritten  using  the  definition  of  G  as 

d  -  6, ,  (G  +  H)  ( 1 2d) 


The  three  basic  equations  which  characterize  the  turbulent  flow  over  the  surface  of  the 
rotating  segment  (the  momentum  integral  equations  (3a)  and  (3b)  and  the  entrainment  equation 
(6b))  are  now  expressed  in  terms  of  the  three  basic  unknown  quantities:  the  streamwise  momen¬ 
tum  thickness  0 ( j ,  the  shape  factor  parameter  H,  and  the  tangent  of  the  angle  between  surface 
streamline  and  external  streamline  directions,  t  =  tan/?. 


NUMTK1C  AI  AN  A I  YSIS 

In  order  to  obtain  a  more  compact  form  of  the  three  basic  equations,  (3a),  (3b),  and  (6b). 
the  expressions  for  the  crosswise  momentum  and  displacement  thicknesses  and  the  streamwise 
displacement  thickness,  equations  (l()a)  through  (l()e),  are  substituted  into  each  of  these  equa¬ 
tions.  After  some  manipulation,  the  single  hyperbolic  equation  in  vector  form  becomes 


aw  aw 

A(W(ln))  —  +  B(W'U.  r?))  —  =  C(W) 

ot,  on 

[lie  unknown  vector  W  is  given  by 
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and  the  known  coefficient  matrices  are 


A(W) 


B(\V) 


C(W)  = 


/ 


•>! 

V:2I 
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m22 

h2 


l2i 

h2 


0n 

3  M  l2 

011 

a,vi,2 

h2 

at 

h: 

3H 

0ii 

d  m22 

011 

3  M22 

)>2 

at 

hi 

3H 

.  hi 

a  d: 

0ii 

a  D-. 

h  2 

if 

h: 

) 


( 1 4h ) 


(Ik  1 


<  1 4d ) 


(  1  ~  —  Ct|  -  Ks0|i  ( I -  M22)  ■  K 1  @  1 1  (2Ms  1 -  O2) 


0, ,  (2Ms,-Ds)  01 


U  h- 


dn 


(  1 4c  I 


e,,(2U)i)  31'  2u>i 

Uh,  04  ‘  "IT 


0nl): 


au 


3 c  1 2 "  20 1 1 X12 1  (777  73  4  K: )  ■  -3  ei i^t 


2enM::  d[_ 

h  2  f  a  n 


(141  | 


K|fl|  |(  IM22  •  D| ) 


1 


-  a  1 1  <  g  •  h  1 )  1 )  ( 


1 

'  I'l 


ai_ 

dl 


K:  I 


•  0ii1-’: 


1  ai_ 

h  2  l  3  n 


K,  ! 


(  1 4  c ) 


14 


Since  equation  (13)  is  hyperbolic,  the  entire  solution  between  asymptotes  could  be  determin¬ 


ed  once  initial  conditions  were  specified  along  these  limiting  lines.  Initial  conditions,  however, 
are  not  readily  available  for  the  geometries  of  interest.  Therefore,  a  solution  procedure  is 
adopted  which  yields  results  only  along  one  side  of  the  bisector  of  the  asymptotes  and  requires 
initial  conditions  along  only  one  asymptote. 

A  first-order  finite  difference  approximation  is  used  to  solve  equation  (13)  in  the  streamline 
coordinates  (4,  p).  The  grid  spacing  is  given  by  A4  in  the  streamwise  direction  and  A p  in  the 
crosswise  direction.  Arclengths  are  approximated  by  the  chordlength;  hence,  reasonably  small  in¬ 
crements  are  required.  Superscript  i  refers  to  the  crosswise  (p)  direction  and  superscript  j  refers 
to  the  streamwise  (4)  direction,  forward  differencing  in  4  yields 

( d  w\  j  v\  j  ♦  ’  u  i 

\  34  /  ~  A4  ( ! 5a) 

and  backward  differencing  in  n  yields 

I d  w\  1  \V»  -  W1  1 

Vah )  ~  ^  ( 1 5  h ) 

If  the  values  of  W ,  { ,  H,  and  tan/3  are  known  at  the  grid  node  ( i , j ) ,  the  values  of  these  three 
variables  can  he  (bund  at  the  node  (i,j  r  I)  by  the  first  order  finite  difference  approximation 


Am 


\vm  *  i  - 


/ \V i.j  -  W'-l.A 


(161 


It  the  A1-1  matrix  is  nonsingular,  equation  (16)  can  be  solved  for  W'-.i  •  I; 


W'-.i  * 

1  -  A4(A'-i)  'em 

-  ^  [( A'-i)  1  B'-i  (W'-J  -  W'd.i)]  4  \yi.j 

Ar] 

d'i 

w  here 

A4 

=  4J  +  1  -4J 

( 1  Sa» 

and 

An 

P1  -  pl_  1 

(IM'. 
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The  streamline  grid  is  shown  on  a  typical  sector  rotating  at  a  constant  angular  speed  t.  in 
Figure  2.  As  experiments  indicate^),  streamlines  (lines  along  which  £  varies)  are  shown  to  ap¬ 
proximate  constant  radius  lines  and  the  crossflow  lines  (lines  along  which  tj  '  aries)  are  shown  to 
approximate  radial  lines.  Streamlines  are  indicated  by  the  variable  j  and  crossflow  lines  by  the 
variable  i.  The  arelengths  and  Arj  are  shown  to  vary  from  point  to  point  on  the  sector  sur- 
•acc.  Grid  nodes,  which  are  defined  as  intersections  of  streamlines  and  crossflow  lines,  are  label¬ 
ed  i,i  for  i  1  to  N,  j  -  I  to  M  where  N  =  number  of  streamlines  and  M  number  ol 
crossflow  lines . 


TRAILING 

EDGE 


STREAMLINES 


LEADING 

EDGE 


CROSSFLOW  LINES 

Figure  2.  Streamline  Grid  on  a  1  epical  Rotating  Sector 
I  igurc  3  shows  the  calculation  sequence  for  the  turbulent  hotmdarv  laver  on  a  t \ pica  1  sec 
:oi.  l  or  convenience,  the  grid  nodes  used  in  the  calculation  procedure  are  numbered  scqucntiallv 
Itom  the  number  1.  It  is  assumed  that  starting  values  of  0|j,  H,  and  t  t an/3  are  given  along  a 
erosstlow  line,  points  numbered  1-7.  Turbulent  boundary-layer  characteristics  are  solved  along 
successive  crossllow  lines  as  indicated  by  the  increasing  numbers  in  Figure  3.  I  or  example,  dat  i 
at  points  |  and  2  are  used  to  obtain  predictions  at  point  K,  data  at  points  2  and  3  ate  used  lot 
predictions  at  point  9,  etc.  As  can  be  seen  from  this  simplified  example  ol  the  calculation  pto 
cedure.  a  fine  grid  is  nccessarv  near  the  hub  of  the  blade  in  order  to  obtain  results  near  the  liatl 
me  edee. 
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A  blade  is  typically  defined  by  a  grid  containing  about  250  nodes.  Experience  to  date  in¬ 
dicates  that  grid  nodes  should  be  distributed  on  the  blade  with  a  greater  concentration  near  the 
inner  radii  and  leading  edge;  and  arclengths  A4  and  Arj  should  increase  uniformly  in  the  direc¬ 
tions  of  the  trailing  edge  and  tip  of  the  blade.  Using  the  Burroughs  7700  high-speed  computer, 
the  boundary -layer  computation  for  this  grid  requires  less  than  10  seconds  computer  tone  ai  a 
cost  of  less  than  one  dollar. 


LEADING 

EDGE 


figure  3.  Points  at  Which  Calculations  Performed 
COMPUTATIONAL  RESULTS  AND  DISC  USSION 

rite  computer  code  for  calculating  the  three-dimensional  turbulent  boundary  layer  lias  been 
evaluated  using  several  geometries  for  which  analytical  and  experimental  data  are  available. 
Among  these  arc  a  flat  plate,  a  two-dimensional  foil,  a  rotating  disk,  and  a  totaling  helical 
blade.  Excellent  agreement  with  test  data  was  found  in  each  case,  flic  turbulent  boundary  layer 
was  computed  for  the  rotating  helical  blade  assuming  both  zero  and  nonzero  crossflow  .  A-  was 
anticipated,  results  of  these  two  computations  showed  agreement  within  15  percent.  Results  arc 
presented  for  only  the  more  complex  of  these  geometries,  the  rotating  helical  blade  with  nonzet 
crossflow.  I  he  turbulent  boundary-layer  development  over  the  surface  of  Blade  1,  a  sector  ol 
300  degrees  included  angle,  was  investigated  experimentally  and  analytically  by 
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Lakshminarayana,  et  al.W.  The  small  pitch  of  the  sector  was  not  included  in  their  analysis,  nor 
is  it  included  in  the  present  investigation. 

The  effect  of  blade  warp  and  skew  is  investigated  using  three  blades:  an  unskewed  Blade 
2A,  a  72  degree  warped  Blade  2B  and  a  72  degree  skewed  Blade  2C.  Blades  2  are  analytically- 
defined  lifting  surfaces  for  which  no  previous  experimental  or  analytical  turbulent  boundary- 
layer  data  exist.  Details  of  the  geometry  and  potential  flow  calculations  are  given  in  Reference 
13.  Initial  conditions  for  this  geometry  are  assumed  to  be  flat  plate  and  rotating  disk  values  for 
turbulent  flow.  The  computed  local  skin  friction  coefficients  are  compared  for  these  blades.  At 
present,  only  the  skin  friction  coefficients  are  computed.  In  the  future,  it  is  anticipated  that  the 
force  and  moment  values  will  also  be  computed,  allowing  a  calculation  of  the  effect  of  viscous 
drag  on  overall  performance. 

Blade  1 

The  turbulent  boundary-layer  characteristics  were  investigated  initially  on  a  single  non¬ 
lifting  rotating  helical  sector  of  92.6  cm  (36.6  in.)  diameter  with  a  300  degree  included  angle, 
hence  of  a  large  chord  length.  Both  experimental  and  analytical  data  are  available  from 
l  akshminarayana,  et  al.(4)  for  comparison. 

1  he  pitch  ol  the  helical  surface  studied  by  l  akshminarayana,  et  ald^l  was  25.4  cm  (10  in¬ 
ches),  lor  a  pitch  diameter  ratio  ot  0.273.  Thus.  the  geometry  can  be  approximated  by  a  flat  cir¬ 
cular  disk  with  a  leading  and  trailing  edge  300  degrees  apart.  The  blade  was  mounted  on  a  22. 86 
cm  (0.75  ft)  radius  hub  and  the  maximum  radius  of  the  blade  was  46.48  cm  (1.525  ft).  The  blade 
was  enclosed  in  a  housing  wall  with  a  0.20  cm  (0.08  in.)  clearance  between  the  blade  tip  and  the 
wall.  The  tree-stream  onset  flow  was  zero  and  the  rotational  speed  of  the  blade,  denoted  to,  re¬ 
mained  constant  at  47  rad  s  (450  rev  min).  The  geometry  of  the  sector  is  shown  in  f  igure  4. 

The  theoretical  analysis  of  I  akshminarayana,  et  al.(t>)  is  based  on  a  rotating  system  describ¬ 
ed  with  cylindrical  coordinates  (r,  6,  /)  where  r  varies  in  the  radial  direction.  6  varies  in  the 
chordwise  direction  ( i . e . ,  circumferential  direction)  and  /  varies  along  the  surface  normal  (axial 
direction).  These  coordinates  are  indicated  in  figure  4  and  are  used  when  comparing  the  ex¬ 
perimental  data  with  the  results  of  the  present  theory. 
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Figure  4.  Geometry  and  C  oordinates  of  Rotating  Helical  Sector 
The  potential-flow  surface  velocity  for  Blade  1  is  simply  ter  along  the  chordlines.  The 
streamlines  are  lines  of  constant  radius  along  which  neither  the  potential  flow  velocity  nor  the 
metric  coefficient  h->  R  varies.  The  momentum  integral  and  entrainment  equations  for  this 
sector  are  thus  simplified  to: 
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the  metric  coefficients  h|  and  h2  are  the  chord  length  and  maximum  radius,  respectively  and  £ 
and  q  are  the  chances  in  the  fraction  of  chord  and  fraction  of  radius,  respectively.  Substitution 


of  these  values  into  equations  (19a)  and  (19b)  yields  the  momentum  integral  equations  given  by 
L.akshminarayana,  et  al.W.  In  their  analytical  study,  only  the  streamwise  momentum  thickness 
0!  [  and  the  limiting  streamline  angle  ft  varied.  A  general  form  of  the  solution  was  assumed  (0|  j 
was  assumed  to  vary  with  radius  and  chord  and  tan/?  was  assumed  to  vary  with  0|  (  only)  and  the 
governing  momentum  integral  equations  solved  numerically  using  the  fourth-order  Runge-kutta 
method.  In  reference  (9),  it  was  assumed  that  the  streamwise  velocity  profile  varies  as  the  1  7 
powet  of  (C  d)  and  a  simple  friction  formula  was  selected.  The  present  method  also  uses  the  cn- 
nninment  equation  and  a  different  solution  method  based  on  assumed  entrainment  and  skin  fric¬ 
tion  relationships  and  assumed  velocity  profiles.  Thus,  some  differences  occur  in  the  predicted 
flow  parameters. 

I  akshminarayana,  et  al.'^l  presented  experimental  and  theoretical  results  at  a  rotational 
speed  of  47  rad  s  (450  rev  min).  The  kinematic  viscosity  v  equaled  1.486  \  l()-s  m-  s  (160  \  I0  (l 
ft--  sec)  which  corresponded  to  air  at  20  C  (68  I  ).  These  conditions,  for  a  Reynolds  number 
based  on  tip  radius  and  tip  rotational  speed  of  7  \  10s  (or  a  0.7-radius  Reynolds  number  of  2.5 
\  106  (see  next  section  for  definitions)),  are  used  in  the  evaluation  of  the  present  method. 

In  ihi.  numerical  solution  procedure,  l.akshminarayana,  ct  al.^l  performed  calculations  ol 
the  laminar  boundary  layer  until  the  flat-plate  critical  Reynolds  number,  Rv  a>r~0  e  4  x 
It'-',  was  reached.  I  he  values  of  the  streamwise  momentum  thickness  0|)  and  limiting  streamline 
angle  /I  vveie  then  used  as  initial  conditions  for  the  turbulent  calculations.  Mow  visualization  ex¬ 
periments  using  a  sublimation  technique  supported  the  use  of  the  flat  plate  critical  Reynolds 
number  to  determine  the  transition  point. 

With  the  present  scheme,  calculations  began  at  the  location  9  0.733  radians  and  \K 

i  R  0.492.  the  grid  node  location  closest  to  the  experimentally-determined  transition  location. 
Initial  conditions  for  the  streamwise  momentum  thickness  fl||  and  the  tangent  of  the  limiting 
streamline  ft  were  obtained  from  the  experimental  data  of  I  akshminaravma,  et  al.*‘b,  but  no 
tadial  variations  were  considered.  Initial  values  were  taken  to  he  0|j  0.24  mm  (0.0008  It). 

tan/J  0.2,  and  H  1.5  at  each  radial  location  along  the  chordvvisc  position  9  0.733 

t  admits. 
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Figures  5  through  7  present  the  results  of  the  turbulent  boundary-layer  calculation  over  the 
surface  of  Blade  ].  Both  the  experimental  and  analytical  results  of  I.akshminarayana,  et  aid9* 
and  Hat  plate  predictions  (adjusted  for  the  region  of  laminar  How  prior  to  8  =  0.733  radians) 
are  indicated  in  the  figures  for  comparison.  The  variation  in  the  stream  wise  momentum 
thickness  On  is  given  in  Figures  5a  through  5c  at  the  radial  locations  xK  =  r/R  -  0.72,  0.82, 
0.93,  respectively.  The  momentum  thickness  is  predicted  to  be  less  than  the  flat-plate  value  given 
by  Prandli \<-7)  approximation 1  (0j|  -  0.036  x  Rv  *'5,  with  U  cut  and  \  r0),  and  less 

than  the  prediction  of  1  akshminarayana,  et  aid9*.  Use  of  the  present  theory  gives  satisfactory 
agreement  with  the  experimental  results. 

The  computed  variation  of  the  limiting  streamline  angle  /)  is  shown  in  Figure  6a  for  r  R 
0.74  and  in  Figure  6b  for  r  R  0.965.  The  peak  in  the  surface  How  angle  at  the  leading  edge  is 
the  result  of  laminar  flow.  The  present  calculations  start  at  9  0.733  radians  (42  degrees),  the 

measured  point  of  transition  at  ;  R  0.72.  These  figures  show  that  the  currently-applied 

analytical  model  predicts  greater  variations  in  the  limiting  streamline  angle  than  observed  for  the 
measured  data .  The  laminar  sublayer  is  not  included  in  the  present  calculations  and  may  be  im¬ 
portant  lot  this  variable.  Values  near  the  tip  may  be  influenced  by  the  annular  wall. 

The  distribution  of  the  local  skin  Irietion  coefficient  for  varying  radial  locations  is  given  in 
Figure  "a  at  9  2.64  radians  and  Figure  7b  at  9  4.76  radians.  Flat  plate  predictions  were 

made  using  the  Sclutenherr*-^1  curve  for  local  skin  friction,  (  j  0.310  |ln-(2Rfj)  >  2ln(2R^)). 
Figure  "  shows  that  the  measured  values  of  the  local  skin  friction  coefficient  are  huger  than  die 
I hit  plate  ptcdicttoMs  and  smaller  than  the  experimental  predictions. 

Blades  2 

In  this  section,  tcstilts  arc  presented  for  the  computed  bouudat  y-lavci  charactciistics  on 
three  analvncallv  -defined  blades*  "P  an  unwarped  propulsor,  Blade  2,\  (similar  to  NNRIX 
Model  43S1).  a  lughlv  waiped  piopulsoi.  Blade  2B  (similai  to  NSRIX'  Model  44‘*S).  and  a 

•  Plaudit's  approximation  was  adjusted  lot  the  region  of  laminar  flow  In  matching  the  values  ol 
laminar  and  turbulent  momentum  thickness  0||  at  the  location  9  0.733  radians. 
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The  geometry  of  the  warped  blade  is  given  in  Table  1.  The  geometry  of  the  unwarped  blade 


is  identical  with  the  exception  that  the  skew  angles  are  zero.  The  blades  rotate  about  a  cylin¬ 


drical  hub  at  a  constant  angular  velocity  cu  in  a  constant  axial  onset  stream  of  speed  V. 


TABLE  1  —  GEOMETRY  OF  ROTATING  BLADE  2B 


Fraction 
of  Radius 

Chord 

Pitch 

Skew  Angle 
(radians) 

Thickness 

Diameter 

Diameter 

Chord 

0.20000 

0.16500 

1 .49500 

0.00000 

0.24000 

0.25000 

0.19700 

1.53500 

0.07854 

0.19800 

0.30000 

0.22900 

1.54400 

0.15708 

0.15610 

0.40000 

0.27500 

1.48700 

0.31416 

0.10680 

0.50000 

0.31200 

1.38900 

0.47124 

0.07680 

0.60000 

0.33700 

1.27700 

0.62832 

0.05650 

0.70000 

0.34700 

1.15600 

0.78540 

0.04210 

0.80000 

0.33400 

1 .04200 

0.94248 

0.03140 

0.90000 

0.28000 

0.93400 

1 .09956 

0.02460 

0.95000 

0.24000 

0.87700 

1.17810 

0.02330 

1.00000 

0.00000 

0.81200 

1.25664 

0.02460 

Fraction 

Thickness 

Chordwise 

of  Chord 

Offset 

Loading 

0.000000 

0.000000 

0.000000 

0.007596 

0.086824 

0.912003 

0.030154 

0.171010 

0.992645 

0.066987 

0.250000 

1 .040899 

0.116978 

0.321394 

1.074103 

0.178606 

0.383022 

1.097916 

0.250000 

0.433013 

1.114882 

0.328990 

0.469846 

1.126317 

0.413176 

0.492404 

1.132938 

0.500000 

0.500000 

1.135109 

0.586824 

0.492202 

1.135109 

0.671010 

0.467219 

1.135109 

0.750000 

0.421900 

1.103578 

0.821349 

0.354824 

0.949248 

0.883022 

0.269637 

0.712632 

0.933013 

0.176665 

0.450324 

0.969846 

0.091629 

0.216717 

0.992404 

0.031655 

0.056755  Number  of 

Blades  5 

1.000000 

1 .010000 

0.000000  Design  Jx 

-  0.889 

% 


Unpublished  experimental  results  at  DTNSRDC  on  a  similar  geometry  indicate  that  the  flow 
over  the  blades  will  be  fully  turbulent  at  a  0.7-radius  Reynolds  number  of  4  \  IOC  A  0.7-radius 


Reynolds  number  is  defined  as: 


iao.7 v V'  *  Pf)2 


(20 
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FRACTION  OF  RADIUS  XD  FRACTION  OF  RADIUS  X 


V 


where:  (c)q  7  =  blade  chord  at  0.7  radius 

V  =  onset  speed 

Jv  =  advance  coefficient,  Jv  =  V/(nD) 

n  =  constant  rotational  speed,  revolutions 

per  unit  time 
l)  -  rotor  diameter 

v  kinematic  viscosity 

The  assumed  flow  conditions  for  turbulent  flow  are  v  -  1.191  x  l()6m-/s  (1.282  \  10'-*  It-  sect. 
V  =  16.9  m/s  (55.3  ft/sec),  and  u>  =  391.5  rad/s. 

No  previous  analytical  or  experimental  boundary-layer  results  exist  for  these  specific  blades 
Initial  conditions  are  assumed  to  be  determined  by  flat-plate  values  for  turbulent  flowl“Jb-4--,|l>, 
with-  0,  |  -  1/  11  0.0142  (Rx)6/7,  C,  =  0.310.  [ln-(2Re)  +  2ln(2Re)j,  and  H  - 
I  [  1  -4.67(C’f) 1  -],  and  by  the  rotating  disk  limiting  streamline  angle  ft  valued),  ft  --  *  11 
degrees  on  the  blade  suction  surface  and  ft  -  -  1 1  degrees  on  the  blade  pressure  surface. 
Because  of  the  uncertainty  in  these  starting  values,  the  computed  results  arc  of  qualitative.  but 
not  necessarily  quantitative,  validity. 

The  computed  potential-flow  streamlines  for  the  blades  arc  shown  in  figure  S.  The 

streamlines  for  the  warped  Blade  2B  and  the  skewed  Blade  2C  coincide  and  both  tire  shown  by 

the  dashed  line  in  figure  8.  Streamlines  at  increased  radial  values  lie  noticeably  inboard  of  the 

constant  radius  lines  on  the  suction  surface.  Except  near  the  tip,  the  same  trend  is  evident  on  the 

pressure  surface.  Streamlines  13  and  14  coincide  on  the  pressure  surface  of  the  three  blades. 

1  he  computed  boundary -layer  characteristics  on  the  three  blades  show  similar  trends  along 

each  streamline,  fherefore.  the  results  will  be  presented  only  along  one  streamline,  streamline 

14.  figures  9  through  II  show  the  computed  distributions  of  0j|.  tan/1,  and  11.  respectively.  In 

creased  values  ot  the  streamwise  momentum  thickness  flj  |  and  the  shape  parameter  II  are 

predicted  for  the  warped  and  skewed  blades  relative  to  the  unskevved  blade. 

figure  12  shows  the  predicted  distribution  of  the  local  skin  friction  coefficient  (  t.  as  o>m 

puled  by  empirical  relationships  given  in  liquations  (8).  lot  Blades  2 A.  2B,  anil  2(  .  flic  sjgmfi 

cant  result  shown  in  this  figure  is  that  the  local  skin  friction  coefficient  is  predicted  to  be  tedne 

ed  on  both  the  suction  and  pressure  surfaces  when  the  blade  is  highly  warped  01  skewed.  Redue 
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tions  of  up  to  30- percent  are  predicted  near  the  trailing  edge.  Similar  reductions  arc  predicted 
for  each  streamline.  Small  variations  in  the  predicted  values  of  the  potential  How  velocity  for  the 
various  blade  designs  do  not  significantly  effect  these  results.  Since  insufficient  data  are  com¬ 
puted  to  permit  a  calculation  of  total  blade  drag,  an  average  reduction  of  15-percent  was 
employed  for  the  average  chordwise  blade-section  drag  relative  to  a  constant  drag  coefficient  of 
0.0085.  Thrust  loading  coefficient  and  power  coefficient  increased  by  one-half-percent  and 
reduced  by  one-percent,  respectively,  hence  an  increase  of  efficiency  of  about  two-percent  for  a 
warped  or  skewed  blade.  The  actual  change  m  efficiency  was  1  percentage  point. 

Several  calculations  were  made  to  identity  the  factors  contributing  to  the  predicted  reduc¬ 
tion  iit  skin  friction  with  increased  blade  warp  or  skew.  Calculations  to  date  seem  to  indicate 
that  these  differences  arise  from  geometrical  effects.  Relative  to  the  unskewed  blade,  the  poten¬ 
tial  flow  velocity  decreases  by  less  than  2  percent  and  the  metric  coefficient  hj  decreases  by  10  to 
20  percent  for  both  the  highly-warped  and  the  highly-skewed  blades.  The  metric  coefficient  Its 
inc; eases  by  Id  tv'  15  percent  for  the  warped  blade  and  by  2  to  3  percent  for  the  skewed  blade. 

A  reduction  is  predicted  in  the  angular  velocity  component  of  b  to  14  percent  for  the  warped 
blade  and  of  only  I  to  3  percent  for  the  skewed  blade,  lhus,  it  does  not  appear  as  though  any 
single  factor  is  responsible  K-,  the  predicted  skin  friction  reduction  with  increased  blade  warp  or 
'kew.  Instead,  the  interaction  of  various  geometrical  factors  seems  to  contribute  to  this 
i  eduction . 

(  <  >\C  I  IT)! NC i  K(  MARKS 

Ibis  paper  presents  analysis  and  results  for  computing  three-dimensional  turbulent  boun¬ 
dary  layers  on  the  surface  of  a  rotating  blade  using  a  momentum-integral  method.  The  analytical 
and  numerical  procedures  for  solving  the  momentum-integral  boundary -layer  equations  have 
been  prevented  and  implemented  in  an  efficient  computer  code.  1  ypical  boundary -layer  com 
putations  require  less  than  It)  seconds  CPU  on  the  high-speed  Burroughs  7700  computer  at  a 
cost  less  than  one  dollar.  The  predicted  boundary-layer  parameters  are  shown  to  give  leasottable 
agreement  with  experimental  data  for  a  simple  body.  Assuming  initial  boundary -layer 
characteristics  determined  by  flat-plate  values  for  turbulent  flow .  increasing  blade  waip  is 


predicted  to  decrease  the  local  skin  friction  coefficient.  1  lie  reduction  in  skin  friction  increases 
the  thrust  coefficient  and  reduces  the  torque  coefficient  yielding  an  increased  efficiency  ol  about 
two  percent  (one  percentage  point).  However,  the  accuracy  of  the  initial  values  and  the  effect  ot 
their  variation  on  the  predicted  results  is  not  known. 

future  work  in  this  area  should  include  an  improved  crossflow  velocity  piofile,  heller 
estimation  of  the  initial  conditions  lor  the  turbulent  boundary -layer  calculation,  and  the  calcul  i 
lion  of  thrust  and  torque  to  predict  drag.  I  \pei internal  boundary -layei  data  arc  needed  to  ) u)l\ 
assess  and  improve  the  present  analytical  work,  future  work  should  also  include  the  prediction 
of  laminar-flow  characteristics  at  the  leading  edge,  laminar  to-iurbulence  iiaiisistiou  pi  edict  ions, 
and  eventually  the  prediction  ol  separation  phenomena  both  on  the  blade  uifaee  and  at  the  to 
where  a  lip  vortex  forms.  Pioneering  work  on  the  tip  legion  has  been  published  by  Shnmioih.  e: 
al.^>> 
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The  turbulent  boundary  layer  results  for  Blades  2  (unskewed  Blade  2A,  72  degree 
warped  Blade  2B,  and  72  degree  skewed  Blade  2C)  presented  in  the  main  text  are  com¬ 
puted  for  a  model-scale  0.7  radius  Reynolds  number  of  4  *  10^.  In  this  appendix, 

similar  results  are  presented  for  these  blades  for  a  full-scale  0.7  radius  Reynolds 

7  -6  2 

number  of  6  *  10  .  The  assumed  flow  conditions  are  V  =  1.191  x  10  m“/s  (1.282  > 

-5  2 

10  ft  /sec),  V  =  16.9  m/s  (55.3  ft/sec),  and  w  =  26.1  rad/sec.  The  propeller 
diameter  is  4.572  m  (15  ft). 

A  finer  computation  grid  was  selected  for  the  full-scale  computations.  More 

grid  nodes  were  added  in  the  radial  direction  in  order  to  compute  the  boundary  layer 

characteristics  closer  to  the  hub  at  the  trailing  edge  of  the  blade.  The  boundary 

layer  results  on  the  model-scale  blades  were  also  recomputed  using  the  finer  grid. 

Both  full-scale  and  model-scale  results  are  given  in  this  appendix. 

Again,  initial  conditions  for  the  streamwise  momentum  thickness  •>  and  the 

sliape  factor  H  are  assumed  to  be  determined  by  flat-plate  values  for  turbulent 
28  29  30 

flow  ’  as  described  in  the  main  text.  The  value  of  the  limiting  streamline 

angle  c  is  also  given  by  the  rotating  disk  value  of  p  =  +11  degrees  on  the  blade 
suction  surface  and  p  =  -11  degrees  on  the  blade  pressure  surface.  Since  no  previ¬ 
ous  analytical  or  experimental  data  exist  for  these  blades  and  the  starting  condi¬ 
tions  are  at  best  a  guess,  the  computed  results  can  be  considered  only  qualitatively 
val id . 


Since  the  local  skin  friction  coefficient  C^.  contributes  to  the  overall  drag, 
it  is,  perhaps,  the  boundary  layer  result  of  most  interest.  Therefore,  the  only 
comparison  presented  is  that  of  the  skin  friction  coefficient  at  the  0.73-radius  for 
tiie  three  blades  at  both  model-scale  and  full-scale.  Figure  13  shows  the  computed 
percentage  change  in  C^.  for  Blades  2B  and  2C  from  the  unskewed  Blade  2A.  The  open 
symbols  represent  the  full-scale  results  and  the  solid  symbols  represent  the  model- 
scale  results.  As  shown  in  Figure  13a,  a  decrease  in  the  skin  friction  coefficient 
is  predicted  on  the  suction  surface  at  both  scales  for  the  warped  and  skewed  blades. 
The  decrease  in  C^.  continues  to  the  trailing  edge  of  the  blades.  There  is  also 
little  variation  in  the  results  for  warp  and  skew  at  each  scale.  Slightly  larger 
changes  in  C^.  are  predicted  for  the  model-scale  blades  than  for  the  full-scale  blade 
Figure  13b  shows  the  changes  in  the  local  skin  friction  coefficient  on  the  pressure 
surface  of  the  blades.  Again,  a  decrease  in  C.  is  predicted  at  both  full-scale  and 
model-scale  for  blade  warp  and  blade  skew.  However,  at  the  trailing  edge'  of  the 


I 


blade,  varies  little  from  the  value  predicted  for  the  unskewed  Blade  2k.  Blade 
skew  is  shown  in  Figure  13b  to  have  a  smaller  effect  on  C^.  than  blade  warp.  Overall, 
Figure  13  shows  that  blade  warp  and  blade  skew  are  predicted  to  decrease  the  local 
skin  friction  coefficient  C^. 


10 

5 

0 

5 

10 

15 

0  0.2  0.4  0.6  0.8  1.0 

FRACTION  OF  CHORD 


- 1 - 1 - 

~1 - ! 

FULL  SCALE 

MODEL  SCALE 

-  °  72  DEG  WARP 

•  72  DEG  WARP  ~ 

0  72  DEG  SKEW 

■  72  DEG  SKEW 

■ 

~  A 

rlR  0.73 

5 

o 

_  i-i 

0 

.  0  0  dn 

• 

.* 

%  •  •iSs-- 

IB 

■ 

1  ! 

_  i  _  i _ i 

(a)  SUCTION  SURFACE 


Figure  13.  F.ffect  of  Warp  and  Skew  on  Predicted  Values  of  Local 
Skin  Friction  Coefficient  C^. 
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